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Abstract 
(N 

^ The phantom force is an apparently repulsive force, which can dominate the atomic contrast of 

a 

^ an AFM image when a tunneling current is present. We described this effect with a simple resistive 

^^ model, in which the tunneling current causes a voltage drop at the sample area underneath the 

^__^ probe tip. Because tunneling is a highly local process, the areal current density is quite high, which 

leads to an appreciable local voltage drop that in turn changes the electrostatic attraction between 
'^ tip and sample. However, Si(lll)-7x7 has a metallic surface-state and it might be proposed that 

Pj electrons should instead propagate along the surface-state, as through a thin metal film on a 

Cd semiconducting surface, before propagating into the bulk. In this article, we investigate the role 

g 
I of the metallic surface-state on the phantom force. First, we show that the phantom force can be 

C observed on H/Si(100), a surface without a metallic surface-state. Furthermore, we investigate the 

1 I influence of the surface-state on our phantom force observations of Si(lll)-7x7 by analyzing the 

T-H influence of the macroscopic tip radius R on the strength of the phantom force, where a noticeable 

> 
OO effect would be expected if the local voltage drop would reach extensions comparable to the tip 

^ radius. We conclude that a metallic surface-state does not suppress the phantom force, but that 

C^ the local resistance Rg has a strong effect on the magnitude of the phantom force. 
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I. INTRODUCTION 

Scanning probe microscopy (SPM) offers the possibility to determine structural and 
electronic properties of a surface on the atomic level [H [2]. The two most common SPM 
techniques are scanning tunneling microscopy (STM) and frequency- modulation atomic 
force microscopy (FM-AFM). With combined STM and FM-AFM, we recently observed a 
so-called phantom force on Si(lll)-7x7 [3J. When the tip is too far from the surface to allow 
the resolution of chemical contrast in the force channel, atomic contrast can still be observed 
in constant-height mode with an applied bias. The resulting images appeared similar to the 
tunneling current images. In our proposed model, the tunneling current is injected from 
the tip into a small area within a radius of ~ 100 pm, the approximate atomic radius of Si. 
This causes an appreciable voltage drop that decreases the electrostatic attraction between 
tip and sample as a function of tunneling current, causing these phantom force images. 
However, a highly localized voltage drop leading to the phantom force effect appears to be 
incompatible with the existence of the metallic surface-state of the Si(lll)-7x7 surface. 
The Si(lll)-7x7 surface is described by the dimer-adatom-stacking-fault (DAS) model [1]. 
One unit cell of the surface consists of 12 adatoms, which have partially filled dangling 
bonds, forming a metallic surface-state [5]. An intriguing question is how electrons 
propagate through the metallic surface-state, which is still under discussion and we refer 
the reader to Refs.[BHT3]. A popular description is that electrons, tunneling from a STM 
tip onto the surface, propagate along the metallic surface-state before entering the bulk. 
To estimate the extension of a voltage drop on a surface with and without a metallic 
surface-state, we used finite element analysis (FEA) software to illustrate these two cases 
[13]. For a surface without a metallic surface-state, we modelled a silicon semiconductor 
sample with a constant conductivity, shown in Fig.flla). To mimic a surface with a metallic 
surface-state, we added a thin metal sheet on top of the sample (Fig.[l]b) -d)). The FEA 
was performed with increasing conductivities of the metal sheet, cr^ = 10^ — , b), cr^ = 10^ — , 
c), and as = 10^ m' '^)' since these metal sheet conductivites cover the range of surface-state 
conductivities noted in Ref . pDJ . with e.g. 10^ m ' 1^ corresponding to 1 ^. In each case, 
a)-d), we defined a highly localized current source on the surface. The current density was 
^^^ ^° (ioo°m)27r ~'^^ i?- ^^ Fig-Pl^)) without a metal sheet, the voltage drop amounts to 
230 mV and is highly focused. In Fig. [lib) - d), with a metal sheet, the voltage drop shows a 
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FIG. 1: (Color online) Finite element analysis of a voltage drop in a plain bulk material a) and a 
bulk material covered with a thin metal sheet on top b) - d) are shown. The bulk material has a 
constant conductivity of at = 10^ (equal to the Si(lll) sample used in the experiment). In b) - d) 
the metallic surface-state of Si(lll)-7x7 is modelled by a metal sheet with a thickness of 100pm. 
The conductivities of the metal sheet increase with as = 10^ — , b), as = 10^ — , c), to 0"^ = 10^ — in 
d). Without the conductive surface layer, a), the voltage drop amounts to 230 mV and is highly 
localized, while the conductive surface layer leads to a reduction and a lateral spreading of the 
voltage drop. 



reduction of its amount and an increasing lateral extension for increasing conductivities of 
the metal sheet. As we have observed the phantom force effect on the Si(lll)-7x7 surface, 
the question remains how the metallic surface-state relates to the phantom force. 
This article gives a description of the phantom force based on our model of an attractive 



electrostatic force in section [IT| Section [ITT] introduces to the equipment and methods 
used for the experiments. Following the conclusions from our finite element analysis, the 
phantom force is expected to occur on a surface without a metallic surface-state. This has 



not yet been demonstrated experimentally. In section IV, we show the phantom force effect 
on a sample that does not have a metallic surface-state. In Section |V] we investigated if 
the metallic surface-state has an effect on our observations on Si(lll)-7x7. If the metallic 
surface-state would establish a constant potential over the whole surface even directly 
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FIG. 2: (Color online) Qualitative distance dependence of the force Fts according to a Morse 
potential, the corresponding force gradient kts and the tunneling current I. For small cantilever 
amplitudes, the frequency shift A/ is proportional to kts- In region I, A/ decreases, whereas / 
increases. In region II, A/ starts to increase. 

beneath the tip, we would expect a delocahzed effect and thus, we would expect the 
observed phantom force to depend on the macroscopic tip radius R, just as the electrostatic 
force between a plate and a semi-spherical tip depends on the radius [IHl HE] . To evaluate 
this, we performed constant height images on Si(lll)-7x7 to extract the ratio between the 
frequency shift due to the phantom force and the tunneling current ('phantom force slope', 
defined in eq. [?]) and relate it to the macroscopic tip radius R, which was determined by 
force versus distance spectroscopy at zero effective bias. 



II. THEORETICAL DESCRIPTION OF THE PHANTOM FORCE 



In this section, we introduce FM-AFM and describe the expected relation between the 
FM-AFM signal and the tunneling current in contrast to the relation between FM-AFM 
signal and tunneling current with a phantom force. Additionally, we mathematically derive 
the contribution of the tunneling current on the electrostatic force. 

In FM-AFM, the forces between tip and sample cause a frequency shift A/ relative to the 
resonance frequency /o of an unperturbated cantilever pTl \TE\ . The cantilever has a stiffness 
k and oscillates with a constant amplitude A at a distance z from the surface. For small 
amplitudes, the relation between the force Fts and A/ is A/ ^ 2^ ' ^ts, where kts = ~^^ 
is the force gradient between tip and sample [19]. For a Morse potential, which describes 



the chemical interaction between tip and sample atom, the force and A/ behave as shown 
in Fig.|2} In region I, as the tip approaches the sample, the force becomes more attractive 
and A/ more negative. Approaching the tip closer to the sample, attractive chemical bonds 
start to form and A/ decreases further. In region II, A/ starts to increase. In this article, 
we focus on region I, which is the region at the onset of current. A more detailed explanation 
of the behaviour between force and A/ is given e.g. in Ref.[l8]. When performing STM, 
the tunneling current exponentially increases with decreasing tip-sample distance [20j. If 
force and current were independent, we would expect, a decrease in the frequency shift as 
the tunneling current increases on a A/ versus / plot. 

However, a surprising characteristic of the phantom force is the increase of the frequency 
shift as the tunneling current increases when plotting A/ against / f3]. 
The phantom force can be modelled by extending the formula of the attractive electrostatic 
force between two metal objects by a tunneling current dependent term. Without the 
influence of the tunneling current we can write 

where V is the potential difference between tip and Cts is the capacity of the tip-sample 
junction. If the tip was a fiat surface A at a distance z to the sample, the derivative of 
capacity with distance would be given by 

The permittivity of vacuum eq ^ 8.85 pF/m can also be expressed as eo ~ 8.85 pN/V^. Thus, 

for A = 20 z'^, a force of about 90 pN would arise for a bias of 1 Volt, increasing with the 

square of voltage. 

The effective bias responsible for the electrostatic force is ^ = Vuas — Vcpd, where Vuas is 

the applied voltage and Vcpd is the contact potential difference between tip and sample [21] . 

While local changes in Vcpd [22] will affect the local electrostatic attraction between tip 

and sample dependent upon V^ias, they cannot explain observations of this phantom force, 

for reasons discussed in Ref. p]: A local change in Vcpd would cause a A/ decrease in one 

bias (assuming the apphed \Vbias\ > \Vcpd\) and an increase in the opposite bias, whereas 

we observe an increase in A/ with significant bias independent of polarity. 

We thus consider the voltage V being modified by a voltage drop caused by the tunneling 



current passing through the sample with resistivity Rs- Therefore, V = Vuas ~ I ■ Rs- The 
electrostatic force is then 

Ft: = -^^(nL - 2VuasiRs + I'Ri) (3) 

with an offset component proportional to V^^^s' ^ term linear with / and a quadratic term 
in /. At typical experimental conditions as in our previous experiments, where V^ias = 1.5 V, 
i?s = 150Mi7 and J = lnA [3], the quadratic term is 5% of the linear term and can be 
neglected (however, for very small tip-sample distances as required for atomically resolved 
AFM on low-conductivity samples, this term can not be neglected). Without the quadratic 
term, equation [3] reduces to 



'' ^ 2 dz 



iy^,,-2VuasIRs). (4) 



In order to determine a relation between the frequency shift A/ and the tunneling current 
/, we first have to calculate the contribution of this electrostatic force to the force gradient, 
kll. After substituting / = Ioe~'^^ into equation and taking the derivative of F// with 
respect to z, equation |4] results in 

1 a Cts -,^2 ( "• ^is dCts 



^r.--^-J^VL,-[^-^^)V,^JR.. (5) 

Since A/ is directly proportional to fcf^, assuming the small amplitude approximation, equa- 
tion [5] can be rewritten as 

A/ = ^"^ VL. -ki'^- "^-) V^AI (6) 



Ak dz^ ""'^ 2k \ dz"^ dz 

The A/ line shows a linear dependence with / with an offset that depends on capacity and 
bias and a slope that is linear with Rs and Vuas- We define, from equation [6} the phantom 
force slope as 

K. Vbias Rs , (7) 



d{Af) /o fd'Cts dCts 



dl 2k \ dz"^ dz 

which is usually expressed in ^. The slope is a measure of the strength of the phantom force. 



III. EXPERIMENTAL METHODS AND SETUP 

The experiments were performed in ultrahigh vaccum (?^ 3-10^^° Torr) and at room 
temperature. The images in this article were all aquired in constant height mode. qPlus 



sensors {k = 1800 — ) were equipped with tungsten tips to probe the sample. The tungsten 

tips were prepared by common techniques hke controlled collision with the sample, field 

emission and explosive delamination |23j . 

The n-doped Si(lOO) samples had a resistivity of 0.008- 0.012 i7cm at 300 K. The surface 

was prepared by repeated cycles of fiashes up to 1250 °C followed by cooling periods in 

the range of several minutes. After cleaning, approximately 300 L deuterium [21] were 

deposited onto the surface at ^450 °C. 

The Si(lll)-7x7 samples used were p-doped with p = 6-9r2cm at 300 K. The surface was 

cleaned by the same fiash routine as described above. 

Investigation of a potential preamplifier artifact 

In the experimental setup, the bias voltage Vbias was applied to the tip with the sample 
referenced to virtual ground via a preamplifier ('preamp'). The preamp is attached outside 
vacuum and amplifies, as a current-to- voltage converter, the / signal by a factor of 10^ ^. 
Since the tip is oscillating, I is an alternating current (AC) with a direct current (DC) offset, 
where only the DC component is measured by the preamp due to its limited bandwidth. 
To investigate if this phantom force effect is not due to a fluctuation of the virtual ground 
of the preamp, we introduced a switch that allows to either connect the sample to real 
ground via a direct ground connection or the virtual ground of the preamp, as schematically 
shown in Fig.jsla). Because the operational amplifler used in the preamp has a limited gain, 
limited bandwidth and a limited slew rate (in contrast to an ideal operational amplifler), 
the virtual ground terminal can deviate from zero, and cross-coupling to the force gradient 
measurement might occur. 

In the upper and lower section of Fig.[3]b), simultaneously recorded / and A/ data are 
presented in constant height mode (switch in position A-B). In the middle section of the 
image, the / signal from the sample was shorted to ground (switch in position A- C). Never- 
theless, the phantom force effect is still present in the A/ signal, which clearly demonstrates 
that the phantom force is not caused by a preamp artifact, but by the current-induced local 
potential deviation outlined above. 
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FIG. 3: a) Schematic of the experimental setup. With an implemented switch, the preamp can 
be used for the / signal (switch in position A-B) or can be shorted to ground (switch in position 
A- C). In b), simultaneous aquired / and A/ data in constant height mode are shown. The images 
were scanned from top to bottom. In the middle section the I signal was shorted to ground to 
check if the phantom force effect is real (i.e. caused by a voltage drop at the sample surface) or 
caused by a preamp artifact. 

IV. THE PHANTOM FORCE ON THE HYDROGENATED SI(IOO) SURFACE 



In this section, we present measurements on the H/Si(100) surface. Exposing Si(lOO) 
to hydrogen saturates the unsaturated dangling bonds [25j. The electronic states of the 
hydrogenated dimers have been shown to be outside the bandgap of bulk Si, meaning that 
in contrast to Si(lll)-7x7, the surface does not have a metallic surface-state [26] . 

In Fig. Ilia), a tunneling current between tip and sample, which leads to the phantom 
force effect, is schematized. Fig.lllb) and c) show simultaneous A/ and / data collected at 
constant height with an applied bias voltage of 1.5 V. The dimer rows can be seen running 
from upper left to lower right. The low contrast is due to our choice of a relatively large 
imaging distance to prevent excessive tunneling currents when scanning over the defect 
area, circled in red. The feature circled in red is most likely a dangling bond, which we 
would expect to observe in A/ data as darker (more attractive). However due to the 
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FIG. 4: (Color online) a) shows a tunneling current, which induces the phantom force effect. 
Simultaneous acquired / and A/ data with atomic contrast are shown in b) and c) for V^ias = 1.5 V. 
In d), the relation between A/ and / data is plotted. Data points with increased / data, as the 
defect outlined in red, show a stronger decrease (less attraction) of the A/ values. However, at 
close tip-sample distances, e), and low bias, in this case 200 mV, the defect appears darker (more 
attractive) , f ) . Image g) shows f ) with low-pass filtering and plane substraction applied for clarity. 
Images are 2nmx2nm, ^ = 100 pm, /c = 1800^, /o = 19131 Hz. 



increase of the tunneling current over it, the phantom force effect causes an increase in A/ 
that makes it appear brighter. To investigate the relationship between A/ and J, we plotted 
in Fig.Hd) the A/ information of each single pixel in image c) versus the corresponding 
pixel of the / information in b). For positive applied bias voltages, the I signal is negative. 



The relation between A/ and / data results in S = -34^, if we assume a linear relation as 

described in equation |6| 

In Fig.llle), the bias voltage is decreased to 200 mV and in order to resolve atomic contrast, 

the tip must be approached to the surface, similar to our previous observations of the 

phantom force [3]. The attractive interaction in the presence of the dangling bond is clearly 

observed in A/ data collected at low bias, as shown in FigJllf). FigJllg) is a low-pass filtered 

and plane substracted image from f) to show the dangling bond with better contrast. 

We demonstrated that the phantom force does not depend on the presence of a metallic 

surface-state and still appears on a sample system as H/Si(100) without a metallic 

surface-state. 



V. THE DEPENDENCE OF THE PHANTOM FORCE ON THE MACROSCOPIC 
TIP RADIUS ON SI(lll)-7x7 

In the following section, we investigate the phantom force on the Si(lll)-7x7 surface. If 
the metallic surface-state plays a role we would expect a delocalized effect. Then we should 
observe a pronounced dependence of the phantom force as a function of the macroscopic 
tip radius R [Ul [IS]. We also discuss the results based on the factors of the phantom force 
slope S introduced in section |TTj 

Figure |5]a) and b) show simultaneously acquired / and A/ data of the Si(lll)-7x7 surface. 
In Fig.jsla) the tunneling current reaches its maximium above the adatoms. In Fig.jslb) the 
brighter adatoms show a repulsive force contribution. The frequency shift is less negative 
over regions with a high tunneling current. A linear dependence between A/ and the / signal 
is shown in Fig.p^c). By fitting the data we extracted a phantom force slope S = 0.73 ^. 

The macroscopic tip can be described by its tip radius R, which we determined by fitting 
the long-range A/ contribution between tip and sample to a model assuming only van der 
Waals interaction [28]. In order to minimize the attractive electrostatic force, we took Af{z) 
spectra while compensating for the Vcpd- Before measuring the A/(Vbias)5 we retracted the 
tip from the sample by 100 pm. This reduced the possibility of tip-sample collisions due to 
drift. The voltage corresponding to the maximum A/ value of the parabolic A/(Vbias) curve 
equals to Vcpd [29]. In Fig|6]a) the Vcpd was determined to 0.4 V. The ^f{z) curves were 
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FIG. 5: (Color online) a) Tunneling current / and b) A/ data taken during constant height 
scanning. A Si(lll)-7x7 unit cell is outlined by a white diamond in a). In b) the adatoms appear 
more bright, which is due to less attraction induced by the tunneling current, c) shows the phantom 
force contribution as a linear dependence between the current and A/ data. The phantom force 
slope E of the line is 0.73^. Images are lOnmxlOnm, A ?a400pm, fc = 1800-, /o = 19130Hz 



and Vf 



bias " 



-1.5V. 



fitted to a model incorporating a parabolic tip shape (as shown in Figj6]b)) in accordance 

to Refs.[2Hl[30]. The fit of the Af{z) data in Figgb) result in i? = 600nm. 

Fig. [7] displays different phantom force slopes S dependent on the respective macroscopic 

tip radius R. Sixteen data points, acquired with two different qPlus sensors (sensor l:red 

triangles and sensor 2: blue dots), are plotted. The data points are widely spread and range 

from radii of 51 nm to 6775 nm. The phantom force slopes vary from 0.51 ^ to 15.74 ^. 

In particular, slopes below 2.0 ^ can be observed for a wide range of macroscopic tip radii 

R. We observe no dependence between the phantom force slope and the macroscopic tip 

radius. This supports our hypothesis of a highly local voltage drop, and suggests that the 

metallic surface-state does not play a role. 

The spread in our measurements of Fig.m (H>5.0^) can be discussed with the aid of 

for d^Cu 

2k^ dz^ dz 



equation g E = -^(^ - ^K)VMasRs- We turn now to the factors -^(^ - ^«) 
and Rs in detail, as Vbias was constant for these measurements. 



2feV dz^ 



The factor —wrA^-h^ — ^^f^) can be calculated, assuming a model of the electrostatic force 
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FIG. 6: (Color online) a) Spectrum of A/(Vbias) for determining the Vcpd between tip and sample 
to 0.4 V. b) Spectrum of ^f{z) taken at Vcpd- The macroscopic tip radius R was extracted by 
fitting the curve to a long range van der Waals force contribution, for b) the extracted R = 600 nm. 
The inset shows the tip radius for a parabolic tip shape. 

Fjf/ of a conical tip (half-angle 6) in front of a metallic surface as described by Hudlet et 
al. [16]. This would be applicable, if the tip and sample surfaces could be modelled by a 
constant potential. Calculations for realistic R= 5 nm and 9= 70°, at conditions summarized 
in Ref. |31]. lead to unrealistic values of S= 68^, much larger than 2.8 ^, the experimentally 
determined average S of the values shown in Fig.[7j We propose that this phantom force 
effect is highly localized. Instead of being described by the macroscopic tip shape, it would 
be better described by a model of the nanoscopic tip cluster. Assuming a plate capacitor 
with C = cq-, we can calculate the phantom force slope using equation 7 and the following 
parameters: /o = 20kHz, k = 1800 N/m, k = 2A~^, Rs = 150 Mfi with an applied bias 
Vbias = —1.5V, at a distance z = AAA and a capacitive area A = (Inm)^. This yields a 
slope S = 2.8Hz/nA, equivalent to the experimental average of 2.8 ^. 
Concerning the factor Rg, we observed in Ref. [3] that the higher the sample resistivity the 
higher the slope H. In our case, the data points with higher H were collected on areas with 
an increased number of defects on the Si(lll)-7x7 surface. The dependence between S and 
the defect density on the Si(lll)-7x7 surface was investigated and is shown in Figjsj For low 
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FIG. 7: (Color online) The phantom force slope H in ^ characterizing the phantom force is plotted 
versus the macroscopic tip radii R. Two different sensors were used, sensor 1 (red triangles) and 
sensor 2 (blue dots) . The phantom force slopes H show no dependence on the macroscopic tip radii 
R. 
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FIG. 8: (Color online) phantom force slopes H plotted versus the defect density on the Si(lll)-7x7 
surface. The phantom force slopes H seem small for less defects on the Si(lll)-7x7 and more 
defects point to increased S. 

defect densities, S seems to be low in contrast to higher defect densities with an increased 
S. But, since the tip shape and Rs changed for each data point, the dependence between 
phantom force slopes and the defect densities is not conclusive and has to be investigated 
in more detail. 

VI. SUMMARY AND OUTLOOK 



We showed in section IV that the phantom force is present on a sample system without 
a metallic surface-state. 
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In section |V} we investigated the influence of a metallic surface-state on the phantom force. 
The experimental observation of the phantom force slope H shows no dependence on the 
macroscopic tip radius R. This infers a highly localized voltage drop and we concluded that 
the metallic surface-state does not play a role for the phantom force effect. 
For a future project we suggest low temperature measurements to investigate the dependence 
of the phantom force on the defect density on Si(lll)-7x 7. In this experiment, the tip would 
be more stable and a controlled exposure of a distinct spot on the surface to e.g. oxygen 
could clarify the dependence between phantom force and sample resistivity. 
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